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Abstract 

We study the flux tube thickness in the confining phase of the (2 + l)d SU(2) Lattice Gauge 
Theory near the deconfining phase transition. Following the Svetitsky-Yaffe conjecture, we map 
the problem to the study of the < ecru > correlation function in the two-dimensional spin model 
with Z2 global symmetry, (i.e. the 2d Ising model) in the high-temperature phase. Using the 
form factor approach we obtain an explicit expression for this function and from it we infer the 
behaviour of the flux density of the original (2 + l)d LGT. Remarkably enough the result we 
obtain for the flux tube thickness agrees (a part from an overall normalization) with the effective 
string prediction for the same quantity. 



1 Introduction 



The distinctive feature of the interquark potential in a confining gauge theory is that the colour flux 
is confined into a thin flux tube, joining the quark-antiquark pair. As it is well known the quantum 
fluctuations of this flux tube (which are assumed to be described by a suitable effective string model) 
lead to a logarithmic increase of the width of the flux tube as a function of the interquark distance 
R. This behaviour was discussed many years ago by Liischer, Miinster and Weisz in [1] and is one 
of the most stringent predictions of the effective string description of confining LGT's. Indeed for a 
non-confining theory one would instead expect a linear increase of the width of the flux tube. 

A natural question is what happens of this picture at the deconfinement point. One would 
naively expect a sudden jump of the flux tube thickness from a log to a linear dependence from the 
interquark distance. However we shall show in this paper that this is a misleading picture. Indeed 
the flux tube width also depends on the finite temperature of the theory. In the standard finite 
temperature setting of LGT's in which the quarks are represented by Polyakov loops this means that 
the flux tube thickness (and its R dependence) also depends on the lattice size in the "time- like" 
compactified direction. 

A tentative answer to this question can be obtained in the effective string framework. As we 
shall see in sect. 2.4 by using a duality transformation it is possible to show that as the temperature 
increases the log behaviour smoothly moves to a linear behaviour, thus excluding a log to linear 
transition at the deconfinement point. However this result strongly relies on the effective string 
approximation (even worse on the gaussian limit of the effective string) and it would be nice to have 
some kind of independent evidence. 

Unfortunately the flux tube thickness (and in particular its dependence on the interquark dis- 
tance) is very difficult to study by Montecarlo simulations. The only existing numerical estimates are 
for the 3d Ising gauge model in which, thanks to the efficiency of the existing Montecarlo algorithms 
for this model, large enough values of the interquark distance could be reached and unambiguous 
signatures of the logarithmic increase of the flux tube thickness (at zero temperature) could be 
observed [2]. 

In this paper we propose an alternative way to address the above question in the vicinity of the 
deconfinement transition using the Svetitsky-Yaffe conjecture which is a very powerful tool to study 
the finite T behaviour of a confining LGT in the vicinity of the deconfinement point, at least for 
those LGT's whose deconfinement transition is of second order. The Svetitsky-Yaffe conjecture [3] 
states that the deconfinement transition of a (d+1) LGT lies in the same universality class of the 
magnetization transition of a d-dimensional spin model with symmetry group given by the center of 
the original gauge group. 

This gives us a non trivial opportunity to check the effective string predictions. If we choose 
a (2+1) dimensional LGT with a gauge group with center Zi (like the gauge Ising model or the 
SU(2) or SP(2) LGT's which all have continuous deconfinement transitions), the target spin model 
is the 2d Ising model in the high temperature symmetric phase for which several exact results are 
known. In particular we shall see that it is possible to study analytically the equivalent of the flux 
tube thickness. Remarkably enough the results that we find agree with the effective string ones 
thus strongly supporting the idea of a smooth transition from a log to a linear behaviour as the 
temperature increases. 

This paper is organized as follows. Sect. 2 is devoted to the discussion of some background material 
on the effective string description of the flux tune thickness and on the Svetitsky Yaffe conjecture. 
In sect. 3 we study the (0-16203) correlator in the 2d Ising model, discuss its asymptotic behaviour 
and extract from it a prediction for the flux tube thickness. Finally in sect. 4 we compare it with the 
effective string prediction and discuss some further features of our analysis. 
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2 Background 



2.1 Definition of the flux tube thickness 

The lattice operator which is commonly used to evaluate the flux density in presence of a pair of 
Polyakov loops (or equivalently in presence of a Wilson loop) is the following correlator 



< P( 0,0)P + (0,R)U p (x ,x 1 , x 2 ) > 
< P(0,0)P+(0,R) > 



< (p{x ,x 1 ,x 2 ,R) >= D ^ n nND+ , n m < U p > (1) 



where P(x\,X2) denotes a Polyakov loop in the spacelike position (2:1,2:2) (in the above equation we 
have chosen for simplicity to locate the Polyakov loops in the positions (0,0) and (0,R) and the xq 
coordinate runs in the compactified timelike direction ) while U p (x) denotes a plaquette located in 
x = (xq, xi, X2). The different possible orientations of the plaquette measure the different components 
of the flux. In the following we shall neglect this dependence (see however the comment at the end 
of the next section). In order to avoid boundary effects we then concentrate on the midpoint of the 
Polyakov loop correlator choosing a generic value of xq and fixing x 2 = R/2. With this choice the 
flux density will be only function of the interquark distance R and of the transverse coordinate x\. 
In the x\ direction the flux density shows a gaussian like shape (see for instance fig. 2 of ref. [2]). The 
width of this gaussian is the quantity which is usually denoted as "flux tube thickness": w(R,N t ) . 
This quantity only depends on the interquark distance R and on the lattice size in the compactified 
timelike direction Nt, i.e. on the inverse temperature of the model (this dependence was implicit in 
the above definition). By tuning Nt we can thus study the flux tube thickness in the vicinity of the 
deconfinement transition 



2.2 Dimensional reduction and the Svetitsky Yaffe conjecture. 

Given a (2+1) LGT with gauge group G and a two dimensional spin model with symmetry group the 
center of G, if both the deconfinement transition of the (2+1) dimensional LGT and the magnetization 
transition of the 2d spin model are continuous then, according to the Svetitsky- Yaffe conjecture [3], 
the two critical points must belong to the same universality class and we can use the spin model as an 
effective theory description for the (2+1) dimensional LGT in the neighbourhood of the deconfinement 
transition. This is the case if we choose for instance SU(2) as gauge group for the (2+1) dimensional 
LGT and Z 2 (the center of SU(2)) as symmetry group of the 2d spin model (i.e. the Ising model). 
Another equivalent choice would be to study the (2+1) dimensional LGT with gauge group Z 2 (the 
Ising gauge model) whose center obviously is again Z 2 . 

In this effective description the Polyakov loops of the LGT are mapped onto the spins of the Ising 
model, the confining phase of the LGT into the high temperature phase of the spin model while the 
combination a(Nt)Nt (where cr(Nt) denotes the finite temperature value of the string tension while 
in the following a with no explicit dependence will denote by default the zero temperature string 
tension) is mapped into the mass scale of the spin model (i.e. the inverse of the correlation length) and 
sets the scale of the deviations from the critical behaviour. In order to describe in the dimensionally 
reduced model the expectation value of eq.Q one must extend this mapping also to the plaquette 
operator. This non-trivial problem was discussed a few years ago by Gliozzi and Provero in [4], 
where they were able to show that in the vicinity of the critical (deconfinement) point the plaquette 
operator of the (2+1) LGT is mapped into a mixture of the energy and identity operators of the 
2d spin model. As a consequence the combination of gauge invariant operators which measures the 
density of chromoelectric flux in a meson is mapped into a three point function < eaa > where in 
the Ising case (i.e. for LGT with a Z 2 center symmetry) = 0. Different components of the flux (i.e. 
different orientations of the plaquette) correspond to different coefficients in the linear combination 
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which relates the plaquette operator of the LGT with the energy and identity operators of the spin 
model. These coefficients will play no role in the discussion and we shall neglect them in the following. 



2.3 Agreement between 2d spin model estimates and effective string predictions. 

If we describe the correlator of two Polyakov loops in the high temperature regime of a confining 
lattice gauge theory in d=3 using the Nambu-Goto effective string we obtain exactly a collection of 
Kq Bessel functions. This was observed for the first time by Liischer and Weisz [5] using a duality 
transformation and then derived in the covariant formalism in [6] . The expression which one obtains 
is: 

oo ( E \ 

<P(0, 0)P(0, R) + ) = Kl 2 2R _!L K h(d _ 3) (E n R). (2) 

n=0 V / 

(see eq. (3.2) of [5]). Notice that in d = 3 the R dependence in front of the Bessel functions cancels 
exactly. 

The argument of the Kq functions is given by the product of the interquark distance R and the 
closed string energy levels: 



E n = aT { l + ^ 



1/2 

(3) 



(see eq. (C5) of [5]). 

It is easy to see that in the large R limit only the lowest state (n = 0) survives and we end up 
with a single Kq function: 

lim (P(0,0)P(0,R) + ) ~ K (E R). (4) 

it— >oo 

where: 

f 1 1//2 

E = aN t \l L . (5) 

Prom this expression we read the Nambu-Goto prediction for the finite temperature dependence 
of the string tension: 

a(N t ) = E /Nt = ah-^j . (6) 

As it is well known this expression cannot be exact since it predicts mean value critical indices for 
the deconfinement transition, however it turns out to be a very good approximation up to rather 
high temperatures (we shall further comment on this below). 

Looking at eq.(@J it is tempting to identify the Polyakov loops correlator in the large R limit with 
the spin-spin correlator in the 2d Ising model which is exactly given by a Kq Bessel function with 
argument mR, m being the mass of the Ising model. This is the origin of the relation mentioned 
above between the mass of the effective Ising model: m and the product a(Nt)Nt of the LGT. Notice 
as a side remark that this correspondence is more general than this particular Ising case, since any 2d 
spin model with a spectrum which starts with an isolated pole has a spin-spin correlator dominated 
at large distance by a Kq Bessel function without prefactors. 

A few comments are in order at this point: 

• The experience with Polyakov loop correlators [7-10] shows that the Nambu-Goto action is a 
good approximation (indeed a very good one) for very large R (much larger than Nt ) and values 
of Nt such that: N t > \jAjo. For higher temperatures (i.e. smaller values of Nt) the deviations 
due to the "mean field " nature of the Nambu-Goto approximation cannot be neglected. 
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It will be useful in the following to define the combination tf = With this definition tf is 

dimensionless and has the meaning of a finite temperature. The Nambu-Goto approximation 
would suggest a deconfinement transition for tf = ~ 1.02 while the bound N t > \J^jo 

mentioned above corresponds to tf < 0.5. It is interesting to observe that the Nambu-Goto 
prediction for the critical temperature is different, but non too far, from the known Montecarlo 
estimates for LGT's: in the gauge Ising model the deconfinement transition occurs at tf ~ 
1.2 [11] while in the ST/ (2) case we have t f ~ 1.13 [12]. 

Looking at Eq we see that it is useful to define the dimensionless coefficient 

(7) 

from which we find 

m = aTp (8) 

The region (tf < 0.5) in which we can trust the Nambu-Goto approximation corresponds to 
the range 0.8 < p < 1 




2.4 Effective string predictions for the flux tube thickness. 

The simultaneous dependence of the flux tube thickness on the two variable R and Nt can be evaluated 
exactly only in the gaussian limit. Including higher order terms in the effective string action makes the 
problem too difficult (even if some recent result in the framework of the covariant quantization suggest 
that some simplification could occur if one chooses to study the whole Nambu-Goto action [6, 13]). 

For the details of the calculations we refer the reader to the paper [2] (see also [14]) For our 
purpose we are only interested in the two asymptotic limits: large Nt and finite R (which is the zero 
temperature limit where we expect a log type behaviour) and the opposite one: large R and small 
Nt which is high temperature limit. 

One finds: 

w 2 ~ -!- log A (N t » R » 0) (9) 

Z7T<7 K c 

As it is easy to see in the second limit the logarithmic dependence is on Nt (the inverse of the 
temperature) and not on R which appears instead in the linear correction. 



3 The 3-point correlators in the 2d Ising model 

As we have seen the study of the width of long color flux tube in the (2+l)d LGT's with gauge group 
Z2, SU(2) or Sp(2) can be translated in the study of the ratio of correlators 

(o-(xi)e(x2)o-(x 3 )) 

(a( Xl )a(x 3 )) { > 

in the high temperature phase of the 2d Ising model in zero magnetic field. 
Since we are interested in the large distance behaviour of such a quantity, we will use the so- 
called Form Factors technique (see [15] for its application in the context of the 2d Ising model without 
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magnetic field). Form factors are defined as suitable matrix elements of an operator <f> between the 
vacuum and an arbitrary n— particle asymptotic state 1 

*&..o»(0i> • • • A) = (010(0)14* [o x ) . . . A an (e n )). (12) 

When the theory is integrable, and the S— matrix is exactly known they can be computed exactly 
as solutions of certain functional equations [16, 17]. In the present case the theory is free and the 
S— matrix is simply S = — 1. 

Such a technique turned out to be an effective tool in order to give approximate expressions for 
the large distance behaviour of the two-point correlators in Integrable QFTs, and recently it has 
been employed to analyze the three-spin correlator in the 3-states Potts model [18]. In brief, it is 
possible to rewrite a generic correlator as a spectral expansion whose building blocks are just the 
form factors. For example, for the connected two-point correlator we have in general 

<0W(O)) C = f; i r M ):f n \F£{eu ...A)| 2 e-"W (13) 

n=1 n - J-oo \ z7r ) 

where E n is the energy of the n— particles state as function of the rapidity. 

Similar expressions can be written for n— points connected correlators, in particular for the three- 
point correlation function we can proceed in close analogy with the analysis of [18]. Hence we can 
write (see figQ): 

/DO JQ JQ 
{F{ f F|(0 12 + io+) e-^cosh^+^cosh^) + 
-DO \ 2 ^) 

1 f°° d6 l d9 2 d93 . *\ \ para n , • ,, n , - h \ .-m(R w cosh B x + R 23 (cosh 2 +cosh0 3 )) 



+ n / Z t. Fx fa (fta) F£(9 32 , 2 x + 31 + ii/>) e" m ^ 13 ™™i + ^<™»"2+™»"3» + 

2 J-00 K 2 ^) 
2 J-00 K 27T ) 

The main ingredients which enter the previous formula are the first few form factors of the operators 
a and e. Let us briefly review their properties (see [15]). 
Spin a and disorder [i operators 

Since we are dealing with a theory of free Majorana fermions, the S— matrix is simply S = — 1. 
As a consequence of the fact that we are in the high-T phase of the theory, symmetry implies that 
the form factors of a and \x are non-zero upon odd and even particle states respectively. 

The VEV of \i is known since a long time and happens to be 

=Fg = B\r\ 1/8 = —-—^m 1/8 = Cm 1/8 , 5 = 1.70852190... (15) 

where we used the exact relation between the coupling constant r (the reduced temperature) and 
the mass of the fermion m = 2ttt. The cluster condition [19] fixes the relative normalization between 
the FF of a and fj,, and implies Ff = (fi). The explicit expression of Fg is given by 

3 g 

FZ(ex,e 2 ,03)=iF? tanh^ (16) 

i<j=i 

and the normalization constant can be fixed by means of the residue equation on the annihilation 
poles 

-iRes ei2=l7r F^(9x,0 2 ,e 3 ) = (1-S)F? = 2F[. (17) 



1 We use the parametrization of energy and momentum in terms of the rapidity 8: E = mcosh#, P = msinh#. 
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Figure 1: Generic configuration for the three-point function (a(xi)e(x2)cr(x3)) . 
Stress-energy tensor © and energy density e 

The FF of the perturbing operator e can be extracted from those of the stress-energy tensor 0. 
Their relationship is given by (we recall that A e = 1/2) 

9 = 4?r(l - A £ )re = 2vrre = me. (18) 

The fact that we have S = — 1 implies that the only non-zero form factor of the trace G is the 
two-particle one F® 

Q 

Ff{6) = -2mm 2 sinh - (19) 

which have been normalized by means of the condition F®(iir) = 2-7rm 2 . The final step is to write 
-F| using the relation Q = me 

1 9 

Fn = —Ff = -2ivrmsinh-. (20) 
m 2 

All the other FFs of e(x) and Q(x) are zero. 
3.1 The width of the flux tube 

We are now in the position to give an analytic estimate for the width of the flux tube in the limit of 
large distances. When the sides of the triangle are large, the leading behaviour of (|14|) is given by 
the first term in the rhs. It is interesting to notice that such a term can be written in an explicit 
way for an arbitrary triangle. The reason lies in the simple form of the 2p form factor of the energy 
operator which allows to write it as follows 

J" ^^ {F ^) 2 FI(0 u + l a + )e~ m ^ cM+R ^ cosh ^ = 

= m(F?) 2 S m?± [°° cosh^i msh ^ e -m(Ri 2 o M + R 23M ) = 

K 11 2 J.^ 2tt 2 2 

= (Ff) 2 sin^±^=^ (21) 



2 VR12R23 
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and hence 

(a( Xl )e(x 2 )a(x 3 )) = (Ff) 2 ~ ^1 = ^ e -™<**+*»> + . . . (22) 

2 -K12-K23 

Since we are interested in analyzing the width of the flux tube at the mid-point between a(x\) and 
a(xz) for large separations of them, we will put R\2 = -R23 = L and i?i3 = 2r in dJ), and we 
define the transverse distance as y = Lcosa+/2 (we also have r = Lsina + /2). Then, we define the 
three-point function in such a geometric configuration as S(r, y) 

S(r,y) = (Ff) 2 ^-^e-W^. (23) 

yi _j_ J.Z 

In order to study the flux tube shape we must study the following ratio 

™ ~ (W (24) 

where (<r(2r)cr(0)) is the two-point correlator between <j(x\) and (j(x 3 ). In the large distance limit 
its leading behaviour is 

/oo in ( zpcr\2 

™( F °f e -2mr coshe = ^L K ^2mr) (25) 
-oo 27r n 

and hence the ratio can be cast in the form 

S(r v) irr e -im^ y 2 +r 2 

P(d,y) = *™ =^—2 k (o \ ■ ^ 
(cr(2r)cr(0)) y z + r z Ko(2mr) 

It is easy to see that for very large separations (mr — > oo) the shape becomes of the gaussian type 
in the transverse variable y 

P(r,y) c -^ Jf^-e-^y\ (27) 

Finally we can study the variance w 2 (r) of P(r,y) wrt y, which exactly corresponds to the width of 
the flux tubes which we were looking for. Let us define: 

/oo 
dyy 2 P(r,y) (28) 
-oo 

as a consequence, setting x = y jr we obtain 

Such an integration cannot be performed exactly, but we can still give its asymptotic estimate in the 
limit mr — > oo. With standard techniques we obtain 

dx ^_ e -2mrVI+^ „ 0Fe- 2mr ( -(mr)' 3 / 2 - — (mr)~ 5 / 2 + 0({mry 7 / 2 )) . (30) 
3 1 + x z \2 32 J 

and combining it with 

g2mr 



K (2mr) 



p 2mr / 1 \ 

— =- I 2(mr) ' + -(mr) " 1/2 + 0((mr)- 3/2 ) (31) 
V 71 " V 8 / 
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we finally obtain 

w 2 (mr) ~ 7T— (l — + Of— -r=)) (32) 

m \ 2mr \{mr) z J J 

which states that the width of the flux tube, in the limit of very large separations between the spins, 
behaves linearly with the separation R = 2r 

w 2 (mR) ~ s + 33 

v ; 2 m 2m 2 v ' 

Finally, it is not difficult to calculate the asymptotic expansion of the higher order momenta 

_ _2n r<x> x 2n 



W W ( r ) = Jll / dx _^_ e -2-r^i+^ (34) 



It turn out to be 



„ (= -2mr- /rN n N N ~ k (-1)3 

w (?n)(r) = ^L e (_L) l im y y l_L r (2£; + j + n + l/2)(mr)- fc - J ', (35) 

v uv y fc=0 j=o 

where the first few orders are as follows 

w ( 2n )(r) = 0Fr(n + 1/2) (-Y (2 + — (n 2 /2 - n - 1/2) + O ( r^—)) . (36) 

\m/ \ mr \{mr) z J J 



4 Results 

4.1 Comparison with the effective string predictions 

Ea. (|33|) is the main result of this paper. If we compare it with the effective string prediction of 
ca. (jlU|) we see a remarkable agreement between the two results. Both show a linear increase with the 
inter quark distance of the flux tube thickness. The linear term appears in both equations with the 
same R/N t dependence and the right dimensions given by i. The factor in front of this correction 
is not the same in the effective string and in the spin model cases, but this is not strange since there 
is a finite renormalization in the mapping between the plaquette operator (in the LGT) and the 
energy operator in the spin model (and similarly between the Polyakov loop and the spin operator). 
Moreover one must recall that the effective string result is obtained in the framework of the purely 
gaussian approximation 



4.2 Isolines of chromoelectric flux 

It is interesting to plot the isolines of the chromoelectric flux as a function of the interquark distance 
in the approximation of the two dimensional Ising model. These can be easily obtained from eq. (|23|) . 
Since in the reduced model we have only one scale m I which combines both the finite temperature 
and the interquark distance of the original LGT these plots can be interpreted in two ways. 

• We may look at them as the result of keeping the temperature of the LGT fixed (i.e. m fixed) 
then as ml increases we are effectively describing an increase of the interquark distance R. It is 
nice to see that indeed as ml increases (from figEJto figEJ) the shape of the flux tube becomes 
more and more narrow, as one should expect for a confining theory. 
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• Alternatively we may think to keep the interquark distance fixed and look to what happens as 
we increase the temperature and approach the deconfinement transition (i.e. as m decreases). 
Again it is nice to see that the flux tube smoothly moves from the narrow shape of a confining 
theory (figEJ) to a shape in which the flux lines become more and more delocalized in the two 
dimensional surface (figEJ- 



4.3 The shape of the flux tube. 

Finally it is interesting to notice that from the explicit expression of the {(J\€2(Jz) correlator one can 
see that the shape of the flux tube is not exactly a gaussian (see ea. (|26j) ). This is indeed also what 
is visible in the montecarlo simulations (see fig. 2 of [2] and the whole discussion at the beginning 
of [20]). Our analysis could suggest a tentative analytic form for these deviations which should 
be valid in the neighborhood of the deconfinement transition and hence extend to this regime the 
functional form valid at low temperature proposed in [20]. This prediction is encoded in the higher 
order momenta reported in eq.s (|35l36j) . It would be nice to check this prediction by performing also 
in the high temperature regime a set of montecarlo simulations similar to those discussed in [20]. 
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Note Added. An interesting application of form factors to the (2+l)-d SU(2) gauge theory [21] 
has appeared at the same time as the present paper. In such a work, the (2+l)-d SU (2) gauge theory 
is generalized to an anisotropic form with two gauge couplings, and the form factors of the currents 
of the SU(2) principal chiral model in (l+l)-d have been used to compute the string tension in the 
anisotropic regime. In this context, the mechanism of dimensional reduction is not related to the 
Svetitsky-Yaffe conjecture. 
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Figure 2: Isolines of chromoelectric flux for m = 1. The circles are centered in correspondence of the 
static quarks and their radius is 1/m. Their extention gives a measure of the region where we expect 
large corrections to the leading behaviour we used to produce both the picture and the estimate of 
the fluxtube width (we recall that we expect that the leading behaviour is valid for mRij >> 1 
where Rij is a generic side of the triangle in figure 1). 
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Figure 3: Isolines of chromoelectric flux for m = 2 
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Figure 4: Isolines of chrome-electric flux for m = 5 




